An Example of Quantum Anomaly in the Physics of Ultra-Cold Gases 
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In this article, we propose an experimental scheme for observation of a quantum anomaly — 
quantum-mechanical symmetry breaking — in a two-dimensional harmonically trapped Bose gas. 
The anomaly manifests itself in a shift of the monopole excitation frequency away from the value 
dictated by the Pitaevskii-Rosch dynamical symmetry [L. P. Pitaevskii and A. Rosch, Phys. Rev. 
A, 55, R853 (1997)]. While the corresponding classical Gross-Pitaevskii equation and the derived 
from it hydrodynamic equations do exhibit this symmetry, it is — as we show in our paper — violated 
under quantization. The resulting frequency shift is of the order of 1% of the carrier, well in reach 
for modern experimental techniques. We propose using the dipole oscillations as a frequency gauge. 

PACS numbers: 67.85.-d, 03.65.Fd 
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Introduction - The current decade is marked by the 
emerging links between uhracold-atom physics on one 
hand and cosmology and high-energy physics on another. 
Examples include kinetics of black holes [lj|, electron- 
positron pair production Q, Zitterbewegunq H, and the 
string theory limits posed on viscosity [j, |5|. How- 
ever, quantum anomalies are usually considered to be 
a purely quantum- field-theoretical phenomenon @, 01 • In 
this article, we suggest a scheme for observing a quan- 
tum anomaly in ultracold two-dimensional harmonically 
trapped Bose gas. 

Quantum anomaly — otherwise known as the quantum 
mechanical symmetry breaking — consists of three ingre- 
dients. First ingredient is an exact symmetry in the clas- 
sical version of the theory in question. Second is a diver- 
gence that appears in the straightforward quantum ver- 
sion of the theory. Third ingredient is a weak violation 
of the original symmetry that emerges in a regularized 
version of the quantum theory. 

The two-dimensional (5-potential has been long recog- 
nized as an example of quantum anomaly in elementary 
quantum mechanics [^, 0]. The classical symmetry of 
the (52-potential originates from the following property: 
under the scaling transformation r ^ Ar, the potential 
transforms in exactly the same way as the kinetic energy 
does. A consequence of this property is the absence of 
any length scale in the corresponding dynamical prob- 
lem, both before and after a straightforward quantiza- 
tion. Next, an analysis of the scattering properties of the 
(52-potential shows a divergence in an all Born 

orders of the scattering amplitude starting from the sec- 
ond. Finally, the subsequent regularization l^, U | leads 
to the appearance of the new length scale a2D- The orig- 
inal symmetry becomes broken, and a quantum anomaly 
emerges. 

Similar anomaly in the case of the potential was 
discussed in Ref. |l2|, along with various regularization 



In Ref. [ISj, Pitaevskii and Rosch predicted a dynam- 
ical symmetry that appears in the classical field the- 
ory (Gross-Pitaevskii equation) of the (52-iiiteracting two- 
dimensional harmonically trapped Bose gas. This sym- 
metry is a direct consequence of the scaling symme- 
try of the 52-potential, described above. The conse- 
quences of this symmetry are (a) absence of the ampli- 
tude dependence of the main frequency (isochronicity) 
and (b) absence of higher overtones in the time depen- 
dence (monochromaticity) of the monopole oscillations 
of the moment of inertia. Both properties were demon- 
strated experimentally [l6| . along with an anomalously 
slow damping, for the case of a very elongated Bose- 
Einstein condensate; its Gross-Pitaevskii equation coin- 
cides with the one for a two-dimensional condensate. 

In this article we address the question of whether the 
Pitaevskii-Rosch symmetry survives quantization. 

In the fully quantized unitary gas, the analogous sym- 
metry has been shown to remain unbroken 17[. The 
quantum correction to the frequency of the monopole ex- 
citation in an elongated condensate [lij was computed in 
Ref. fisl 

The zero-temperature equation of state of the 
two-dimensional Bose gas.- The low-density zero- 
temperature quantum field theory (QFT) expression for 
the chemical potential of the two-dimensional Bose gas 
is well known [T^, [2^ : it reads 
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where 020 is the two-dimensional scattering length, 
7 = 0.5772 ... is the Euler's constant, n is the two- 
dimensional density, and x{z) = _^ \(-z) ■ ^'^^ ^ given 
two-dimensional two-body interaction potential V(r), 
its scattering length 020 is defined as the radius of a 
hard disk whose zero-energy s-wave scattering ampli- 
tude equals the one for the potential V. Here, W-i{z) 
is the minus- first branch of the Lambert's W- function 
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2l| . For the small values of the gas parameter na^D, 
the factor % is a logarithmically slow function of the 

density: l/ln(l/z) + Cl(ln(ln(l/0))/ln(l/z)2). 

The expression ([T]) is an inverse of a more conventional 
formula n{fi) — {m^,/ Anh^)\n{A1i^ / e'^'^^^miia\jj), shown 
to be the leading term in an expansion in powers of 
l/XniAh^ /e^'^+^ racial jj) 

In the case of three-dimensional short-range- 
interacting atoms tightly confined to a two-dimensional 
plane by a one-dimensional harmonic potential, the 
two-dimensional scattering length can be expressed 
through the three-dimensional scattering length 
azD and the confinement size = -\/2?i/ (mujz) 

as a2D = C'2Da^xpJ-(v^/2)(az/a3D)], where 
C2D = 1.4795... 0, [Ij- The chemical poten- 

QFT 



tial becomes fJ,{n) 
g2D = ATTth^ /m is 
coupling constant, e = 
rameter governing the 



g2Dn{xicr{n)e ^)/e) where 
the "bare" two-dimensional 
a-ZD I \fT^o-z is the small pa- 
proximity to the classical 
In the limit 

e ^ min (1, 1/1 ln((T(n))|), the factor x approaches 
the density- independent constant e, and the chemical 
potential converges to the prediction of the classical field 
theory (CFT), otherwise known as the Gross-Pitaevskii 
equation; there, the chemical potential reads 



limit, and (j{n) — ne'^^'^^{C2D)'^nd1. 



^{n) = g2D n 



(2) 



Hydrodynamic equations - The zero-temperature hy- 
drodynamic (HD) equations for a two-dimensional har- 
monically trapped Bose gas read 



(3) 



g-^v + {v- V,)v = -{l/m)VMn) + ^Ho(r)] , (4) 



where m is the atomic mass, n = n(r, t) is the atomic 
density, v = v{r, t) is the atomic velocity, T4jo(r) = 
is the trapping potential energy per particle, 
is the chemical potential, and r — xe^ + ySy. As- 
suming that the gas is a Bose condensate with no vor- 
tices, the velocity field can be assumed to be irrotational: 
V = ?iV$/m, where $(r) is the potential of the velocity 
field. Under this assumption, the HD equations pi4p 
can be rewritten in a Hamiltonian form, ^n{r, t) — 

j-[H,nir,t)]uB , i*(r, i) - j-[H , , t)]^^ ■ The 
Hamiltonian is represented by a sum of two parts, the 
first being in turn a sum of the kinetic and interac- 
tion energies and the second being the trapping energy: 

H = Ho + Hw neve Ho = Jd\ [^n{V^f + ein)^ 
Huo = /d'^r VHo(r), the "hydrodynamic commutator" 
[•, 'Ihd is given by the Poisson brackets with respect to 



the (n, <I>) canonical pair, 

[A(n(r, i),$(r, t)), B{n{r' , t),$(r', i))]HD = - J x 
|-^(n(r, t),$(r, t))^ (^^W^, t)M^, t)) 
^A{n{v, t)M^, t))^ (I^^W^' ^)'*(^' ^)) 

and e{n) = f^dn' /i(n') is the microscopic energy density. 

The Pitaevskii-Rosch symmetry and the quantum 
anomaly at the HD level- Introducing the generator of 
the scaling transformations, Q = ^ /(i^rn(r • V$) (see 
Ref. [31 for example), one obtains the following set of 
commutation relations: 



QFT d 

[Q, i?o]HD = 2iHo + ia2D- Hq 



da2D 



(5) 



[Q, ^hoIhd '^^^ -2zffHO, [i?HO, i^olHD ^= ^t^'Q. No- 
tice that the classical field theory chemical potential ([2]) 
does not depend on the two-dimensional scattering length 
a2D- Then, the commutator ([S]) becomes 



QFT . 



[Q, Hohu ""^^ 2iHo 



(6) 



In this case, the observables Hq, Huo, and Q form a 
closed three-dimensional algebra, identical to the one 
discovered by Pitaevskii and Rosch [15] at the Gross- 
Pitaevskii equation level. 

However, the more accurate quantum field theory pre- 
diction for the chemical potential ([T]) depends explicitly 
on a2D- According to the Eqn. [5l the classical commu- 
tation relation ^ becomes corrected by ia2D gf^^ Hq. 
Since the correction term is not generally expected to 
be a function of three original members of the algebra, 
the algebra opens; such an opening constitutes a quantum 
anomaly. 

Castin-Dum-Kagan-Surkov-Shlyapnikov equations.- If 
the factor x in the chemical potential expression ^ was a 
constant, the HD equations pi4p could be solved via the 
Castin-Dum-Kagan-Surkov-Shlyapnikov (CDKSS) scal- 
ing ansatz [1^ [2^. Note however, that x is a very slow 
function of the density, and thus the scaling ansatz should 
approximately hold. Consider the following ansatz: 

"(^' t) = i^no[ l - {r/iX{t)R)n Hr, t) - e^r^, 
where Rtf = \/'2g2D{no)no/mu!'^ is the steady-state 
HD curvature of the density distribution (that also cor- 
responds to the Thomas-Fermi radius of a gas with a 
density-independent coupling constant fixed to 521) ("reo)), 
g2D{n) = (dfj,/dn) is an effective density-dependent cou- 
pling constant, no is the steady-state peak density; the 
scaling parameter A obeys the CDKSS equation 



A = 



A3 



(7) 
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with u(A) = g2D{no / X^) / §20(^0) ■ It can be shown that 



the above ansatz solves the HD equations (I3l4p almost 
everywhere, with the exception of an exponentially nar- 
row ring close to the edge of the cloud. 

For example, let us define the thickness AR{S) of 
a ring of positions r such that the effective coupling 
g2D{n{r)) differs from that in the center by a relative 
correction 5 or more, but the density is still greater 
than zero: \g2D{'n{r)) / g2D{nQ) — 1\ > 5 and r < Rtf- 
(Here, Rtf is the "true" Thomas- Fermi radius defined 
by T4jo(-Rtf) = M(^io)-) Then one can show that in 
the limit e — > 0, (j{n) — const, the leading behavior of 
the thickness Ai?(5) is Ai?(5) = (1/2) exp(-(5/e) i?Ti.. 
Likewise, one can show that the curvature radius Rtf 
is exponentially close to the Thomas-Fermi radius of the 
cloud: Rtf ~ Rtf- 

Anomalous frequency shift of the monopole frequency- 
Linearization of the equation ([T]) for small excitation am- 
plitudes readily gives the frequency of small oscillations, 
n = rJo(l + 5(ARa=0)), where ^ 



,5(ARa = 0)) = 



1 x(cr(no)e 
4 (l-x(a(no)e-7))2 



(8) 



The deviation of the monopole frequency f2 
from the classical field theory prediction Hq = 

2lj is a manifestation of 



e<min(l,l/|ln(CT(n))|) 



the quantum anomaly. 
Here and below, AR^ 

the monopole oscillations. 
For e < min | 1, 



(A„ 



)/(A max \ 



|ln(a(«o))| 

correction converges to 
J(ARa=0)) 



1) is the aspect ratio for 



the relative anomalous 



1 flSD 



(9) 



In particular, this estimate shows that the effect of the 
anomaly can be enhanced using a Feshbach resonance. 

Furthermore, a numerical analysis shows that un- 
der quantization, the monopole frequency becomes 
amplitude-dependent. Fig. [1] shows the corresponding 
prediction of the CDKSS equation ([7]) compared to both 
the results of the full HD treatment p|4p and the limiting 
values 

A possible experimental scheme for detecting the 
anomalous frequency shift- Consider the following ex- 
citation scheme. Initially, the cloud is prepared in the 
ground state of a frequency Winit. trap at a position 
X = y — Q. Att = 0, the trap is simultaneously 
relaxed to a lower frequency u and shifted to a new 
position X — xq, y = 0. This initial condition will 
induce a superposition of a monopole oscillation that 
starts from the lowest in cloud size turning point and 
a dipole oscillation — whose frequency ui does not de- 
pend on either interaction strength or the oscillation 
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FIG. 1: (color online). Relative anomalous shift of the 
monopole frequency, 5(ARa) = {fl — f2o)/f^o as a function 
of the amplitude of the excitation, for different values of the 
transverse confinement frequency ljz- Solid lines: predic- 
tion of the modified Castin-Dum-Kagan-Surkov-Shlyapnikov 
(CDKSS) model ((T)). Crosses: prediction of the HD equa- 
tions (pitt . Arrows: the a^o < limit dSj) of the CDKSS 
model prediction ^ for the zero-amplitude shift. The rest of 
the parameters corresponds to A*' = 10^ rubidium 87 atoms 
{asD = 5.2 nm) confined longitudinally in a frequency uj har- 
monic trap, where u = 2tt x 41 Hz for the uj^ — 2-k x 45 kHz 
data and uj = 2n x 40 Hz for the rest of the data. The steady- 
state peak density no entering the CDKSS model was derived 
from the number of atoms TV using the Euler equation Q. 
An possible experimental setup corresponding to the 45 kHz 
curve is described in detail in the concluding section of this 
article. Insert shows the monopole vs. dipole beat signal for 
the optimal choice of the initial trap shift (see main text). 
Parameters correspond to the 45 kHz curve taken at the as- 
pect ratio ARx — 0.486. We also show the initial linear part of 
the slow envelope, (A^„^/6 + ((A^„^ - l)/(12A^<j^)) 5(ARa = 
0.486) Q,ot plotted with a ±10% uncertainty in the value of the 
anomalous correction (5(ARa) = 0.00760. The CDKSS model 
was used for calculations. 



amplitude — that starts from the left(right) turning point 
for a;o > 0{xo < 0). 

What we suggest is to measure the spatial mean of 
the square of the horizontal displacement with respect to 
the new trap center, s{t) = J dxdy{x — xo)^n{x, y, t). 
After a lengthy but straightforward calculation it can 
be shown that this observable will evolve in time as 
s[t)/R^p = P + Acos[17ot] + Scos[rJo(l + '5(ARA))i + 0], 
where P = (Ai„., + 1)/12A2„„, + A = 77V2, 

B = {Xtax - l)/12A^<j^, = TT, 77 = xo/Rtf, and 
Rtf is the effective Thomas- Fermi radius for the fi- 
nal, frequency lu trap (see a definition above). Choos- 

B. It is easy to 



mg 1] 



leads to A 



see that in this case, at t — the beats between the 
dipole and monopole oscillations have a node, with the 
first crest reached in 1/{2S{AR\)) monopole periods. In 
terms of the initial trap, the optimal shift of the trap can 
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be expressed as = C^™'*- ) — 1 ^tf, init., where 

Rtf, init. — v^w/winit. i?TF is the Thomas-Fermi radius 
for the initial, frequency Winit. trap. Note as weh that 
the right turning point for the oscillations of the scaling 
parameter can be expressed as Xmax — ^/winit./w. In our 
excitation scheme, the monopole oscillations start from 
the left turning point that in turn corresponds to the low- 
est cloud size. Insert of Fig. [1] presents an example of the 
projected beat signal for a typical set of parameters. 

The present in any realistic trap anharmonicity and 
anisotropy may complicate the matters. For the an- 
harmonicity, one can show that for a quartic correc- 
tion of a form (1 -|- Q;r^/2ag) 
(flo — ^yh/muj being the size of the ground state), the rel- 
ative (to the monopole frequency) shift of the monopole 
frequency is (5anharm. = jfaV Ne, to the leading or- 
der in both a and the amplitude. The doubled dipole 
frequency — that serves as a reference — is unshifted in 
that order. On the other hand, the anisotropy leads to 
a splitting of the dipole frequency with the monopole 
frequency situated in between the two resulting frequen- 
cies. If one chooses the doubled X-dipole frequency 2uJx 
as a reference, the relative to this frequency monopole 
shift will be given, to the lowest order in both anisotropy 
uJx-u}y and the amplitude, by (5anisotr. = ii^x-^y)/'2ujx,y, 
where uj^^y is any of the two frequencies. Both the anhar- 
monicity ((5anharm.) and anisotropy ((Janisotr.) shifts must 
be kept below 1% to allow for observation of the quantum 
anomaly. 

Vortex-antivortex pair creation has been suggested 
as the dominant mechanism for damping of the two- 
dimensional monopole oscillations Q, all the conven- 
tional channels being suppressed due to the Pitaevskii- 
Rosch symmetry. Here, the amplitude of the oscilla- 
tions decays as \max{t) = Xmax{t)/{1 +t/Ty/^°, where 

_ _ 737280 /^ "init. N^ { A„i„(0) y° ln(4mg2D jV/yTFri'') „ 
^- 127 I 1^ ; \\^aAO) J {mg2D/h^)'>^*VN ■ ^' 

it is assumed that the amplitude of oscillations is large 
(i.e. Xmax ^ 1) and the gas is dilute (i.e. <C a^)- 

Let us now propose a concrete example of a trap suited 
for an observation of the quantum anomaly. The num- 
bers proposed are those of rubidium 87 in its 5Si/2,F = 
2,mF = 2 ground state. A blue detuned standing wave 
(wavelength 532 nm, power P, waist Wq) crosses, at a 
right angle, the symmetry axis z of a quadrupolar mag- 
netic field; the magnetic field gradients are given by b' in 
the xy plane and —2b' along z. This results in a series 
of pan-cake traps: there, the magnetic field is responsi- 
ble for the weak horizontal trapping, while the standing 
wave gives a strong vertical confinement. Atoms could 
be loaded in one of the nodes of the standing wave, at 
a distance z — d from the trap center; there the mag- 
netic field is bounded from below by Bq = 2b'd. By 
construction, this trap should be isotropic in the horizon- 
tal plane. Furthermore, the residual anisotropy defects 



could be further compensated by adding small magnetic 
gradients in the horizontal plane. For b' = 167G/cm, 
d = 400 /im, P = A W, and wi^ = 500 /im, the oscil- 
lation frequencies are 41 Hz x 41 Hz x 45 kHz, and the 
magnetic field minimum Bq — 13.4 G well prevents any 
Majorana losses. The photon scattering rate is as low as 
2.6 X 10"'^ s~^, corresponding to 30,000 monopole mode 
oscillations. For N — IQ'^ atoms (ten times the number 
used in Fig. [T]), the anharmonic parameter is bounded 
from above by a = —2.1 x 10~^. The predicted Thomas- 
Fermi radius of 15 //m is much smaller than both d and 
Wo; accordingly, the anharmonic shift has a negligibly 
small value of Janharm. = 5 X 10~^. The chemical po- 
tential oi fx/h = 1.7 kHz is well below the transverse fre- 
quency; this ensures that the trap is well in the two- 
dimensional regime. Finally, taking oscillation aspect 
ratio of ARa — 0.486 as an example, one obtains the 
damping time r which is 1576 times longer than the pe- 
riod of the monopole oscillations. Overall, these figures 
bring our proposal within reach of modern experimental 
technology. 

Summary and outlook- In this article, we propose 
an experimental scheme for observation of a quantum 
anomaly in a two-dimensional harmonically trapped Bose 
gas. The effect consists of a shift of the monopole exci- 
tation frequency away from the value dictated by the 
Pitaevskii-Rosch dynamical symmetry [15i] . The shift we 
predict is only of the order of 1% of the base frequency. 
To detect it, we propose using the dipole oscillations as 
a reference frequency. 

We are grateful to Eric Cornell, Steven Jackson, and 
Felix Werner for enlightening discussions on the subject. 
Laboratoire de physique des lasers is UMR 7538 of CNRS 
and Paris 13 University. LPL is member of the Institut 
Francilien de Recherche sur les Atomes Froids (IFRAF). 
This work was supported by grants from the Office of 
Naval Research iN00014-06-l-0455) and the National 
Science Foundation {PHY-0621703 and PHY-0754942). 



[1] L. J. Garay, J. R. Anglin, J. I. Cirac, and P. ZoUer, Phys. 

Rev. Lett. 85, 4643 (2000). 
[2] P. O. Fedichev, U. R. Fischer, and A. Recati, Phys. Rev. 

A 68, 011602 (2003). 
[3] J. Y. Vaishnav and C. W. Clark, Phys. Rev. Lett. 100, 

153002 (2008). 

[4] C. V. Johnson and P. Steinberg, Physics Today 63, 29 
(2010). 

[5] J. E. Thomas, Physics Today 63, 34 (2010). 

[6] S. B. Treiman, R. Jackiw, B. Zumino, and E. Witten, 
Current Algebras and Anomalies (World Scientific, Sin- 
gapore, 1985). 

[7] J. F. Donoghue, E. Golowich, and B. R. Holstein, 
Dynamics of the Standard Model (Cambridge Univer- 
sity Press, Cambridge, U.K., 1992). 

[8] R. Jackiw, in M. A. B. Beg Memorial Volume , edited by 



5 



A. All and P. Hoodbhoy (World Scientific, New York, 
1991). 

[9] B. R. Holstein, American Journal of Physics 61, 142 
(1993). 

[10] L. R. Mead and J. Godines, American Journal of Physics 
59, 935 (1991). 

[11] S.-L. Nyeo, American Journal of Physics 68, 571 (2000). 
[12] H. E. Camblong, L. N. Epele, H. Fanchiotti, and C. A. 

Garcia Canal, Phys. Rev. Lett. 87, 220402 (2001). 
[13] H. E. Camblong, L. N. Epele, H. Fanchiotti, and C. A. 

Garcia Canal, Phys. Rev. Lett. 85, 1590 (2000). 
[14] S. A. Coon and B. R. Holstein, American Journal of 

Physics 70, 513 (2002). 
[15] L. P. Pitaevskii and A. Rosch, Phys. Rev. A 55, R853 

(1997) . 

[16] F. Chevy, V. Bretin, P. Rosenbusch, K. W. Madison, and 

J. Dalibard, Phys. Rev. Lett. 88, 250402 (2002). 
[17] Y. Castin, Comptes Rendus Physique 5, 407 (2004). 
[18] L. Pitaevsku and S. Stringari, Phys. Rev. Lett. 81, 4541 

(1998) . 



[19] V. N. Popov, Functional Integrals in Quantum Field 
Theory and Statistical Physics (Reidel, Dordrecht, 1983). 

[20] C. Mora and Y. Castin, Phys. Rev. A 67, 053615 (2003). 

[21] W--l{z) is the real W < -1/e solution oi z = We^ in 
the range — 1 < 2: < 0. 

[22] C. Mora and Y. Castin, Phys. Rev. Lett. 102, 180404 
(2009). 

[23] D. S. Petrov and G. V. Shlyapnikov, Phys. Rev. A 64, 
012706 (2001), note that C2D = 1.4713 ... in this article 
is slightly different from C2D ~ 1.4795 . . . obtained in 

[24] L. Pricoupenko and M. Olshanii, J. Phys. B 40, 2065 
(2007). 

[25] Y. Castin and R. Dum, Phys. Rev. Lett. 77, 5315 (1996). 
[26] Y. Kagan, E. L. Surkov, and G. V. Shlyapnikov, Phys. 

Rev. A 54, R1753 (1996). 
[27] Note that we have used the property z-j^xi^) = 

X{z)'^/{1 — x(^)) that follows from the properties of the 

Lambert's W-function. 



